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An extension of a method of Chen and Lahlum 1 for estimating the 
distributions of high rain rates is described. Through application of the 
statistical theory of extremes to available yearly maximum rain-rate data, 
a reasonably accurate distribution is obtained. The calculated results 
agree well with previously obtained 20-year data. Application of this 
method to Weather Bureau Rainfall Intensity-Duration-Frequency Curves 
has yielded 50-year distributions of 5-minute rain rates for 36 locations 
in the western United States; others can be similarly obtained. Such 
long-term rain-rate distributions are valuable for microwave radio-path 
engineering, especially in the western United States where shorter-term 
data sources are inadequate. 

I. INTRODUCTION 

References 2 and 3 describe a procedure for obtaining long-term 
( ^ 20 years) distributions of 5-minute rain rates from data published 
by the National Climatic Center. 4 Such distributions have been ob- 
tained for 202 locations in the eastern and midwestern United States 
and applied to path engineering of 11-GHz radio links. 

However, the excessive short duration rainfall data 4 on which these 
distributions are based contain only rainfalls that exceed an excessive 
rainfall threshold denned by the National Climatic Center. 23 For 
example, the threshold is 75 mm/hr for 5-minute intervals. In low 
rain-rate areas, such as Oregon and Washington, almost all rainfalls 
do not exceed the threshold and, hence, are not included in the excessive 
short-duration rainfall data. For example, at Spokane, Washington, 
the 5-minute rain rate exceeded the 75 mm/hr threshold only once in 
the 20-year period from 1953 to 1972. This data source, therefore, is 
an unsatisfactory basis for radio-path engineering in such areas. On 
the other hand, processing other longer-term data — say 50 years — is 
tedious and costly. This has motivated the search for an alternative 
method. Fortunately, the statistical behavior of the extremes of a 
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random variable has been extensively investigated. 6-10 This paper de- 
scribes a method for obtaining distributions of high rain rates by 
applying this theory to the yearly maximum rain-rate data published 
by the National Climatic Center. 

In an unpublished work, Chen and Lahlum 1 have applied the 
theoretical distribution of yearly maximum 5-minute rain rates and 
an empirical extrapolation to obtain the rain-rate distribution in the 
range of interest to radio-path engineering. In this paper, we extend 
Chen and Lahlum's method by incorporating the theoretical distri- 
butions of yearly kth largest 5-minute rain rates for A; ranging from 1 
to 12. The application of the higher-order statistics of extremes elimi- 
nates the need for empirical extrapolation. 

In this paper, a 5-minute rain rate corresponds to the average value 
of the randomly varying rain rate in a 5-minute interval and is calcu- 
lated as AH/t, where AH is the 5-minute accumulated depth of rainfall 
and r = 1/12 hour = 5 minutes. For illustration, only the statistics of 
5-minute rain rates are discussed in this report. The method is also 
applicable to other integration times. 
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Fig. 1 — Comparison of 20-year rain-rate distribution calculated by extreme 
statistics method (dashed line) with 20- and 22-year data (solid lines) for Newark, 
New Jersey. The difference between the 20- and 22-year data also indicates the 
instability of high rain-rate statistics. 
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Fig. 2 — Comparison of 20-year rain-rate distribution calculated by extreme 
statistics method (dashed line) with 20- and 22-year data (solid lines) for Wilmington, 
North Carolina. The difference between the 20- and 22-year data also indicates the 
instability of high rain-rate statistics. 

II. THE STATISTICS OF EXTREMES 

Let R be the randomly varying 5-minute rain rate and Rk be the 
fcth largest 5-minute rain rate in a given year. In other words, R\ is 
the yearly maximum 5-minute rain rate, R2 is the yearly second largest 
5-minute rain rate, etc. The value of Rk varies from year to year, and 
the probability distribution of Rk is the subject of the statistics of the 
fcth extreme. 

Many sets of rain-rate data 11-15 indicate that rain-rate distributions 
in the moderate and low rain-rate region can be closely approximated 
by the lognormal distribution. In the tail region, the time bases are 
usually insufficient to yield stable results for testing the lognormal 
hypothesis. Figures 1 and 2 show instability at extreme values occurring 
in 20-year time bases. We will assume that the rain-rate distributions 
are lognormal and proceed to show that the calculated distributions of 
extreme rain rates agree well with the data as displayed in Figs. 1 
through 9. 

Let 

x = In R, (1) 



RAIN RATE 1113 



CO CD 
in < 



10'- 



10' = 



co < _ 



10 u - 



- 




1 


LA GUARDIA AIRPORT. NEW YORK 


■ 


20 YEARS (1953-1972) 






= S-12— \ 




- S = 8--^ 




- S = 4 


"^ V 


— S = 2 — 


""" ""*• ^\ 


— 


X \ 


- 3-1- 


--^\ 




N\\ 


_ 


^ V \ 




V> ^v 


— 


^v 




^\ 




N ^V 


— 


"*N. 


~ 


*\ 


- 


vv 








\\ 




^x 




>^ , 


1 I 


I I I 1 I I 1 I I Nl I I 



20 40 60 



80 100 120 140 160 180 200 220 240 260 280 300 
5-MINUTE RAIN RATE IN MM/HR 



Fig. 3 — Comparison of 20-year rain-rate distribution calculated by extreme 
statistics method (dashed lines) with the 20-year data (solid line) at La Guardia 
Airport, New York, New York. 



and 



x k = In R k . 



(2) 



The lognormal hypothesis for R is equivalent to the statement that x 
is (approximately) normal. The distribution of the fcth extreme, Xk, 
as derived by Cramer 8 is 



where 



P(xk ^ X) = 1 - e-'~ v T, 

N=0 



y = a(X- U) 



k-l P -Nv 



N\ 



(3)* 



(4) 



and is called the reduced variate. In this expression, a and U are scale 
and location parameters, respectively, and are related to the sample 
mean and sample standard deviation of X\. Notice that the distribution, 
P(xk ^ X), for any k, is completely determined by the two parameters 
a and U. These two parameters can be calculated from the measured 



* The cumulative distribution function (3) is obtained by integrating the prob- 
ability-density function derived by Cramer in Ref. 8. 
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yearly maximum rain-rate data. Let 

R i (j),j= 1,2,3, ...,Jf 



(5) 



be the measured yearly maximum 5-minute rain rate in M years of 
measurements. For example, Ri(7) represents the yearly maximum 
5-minute rain rate observed in the seventh year of an M year experi- 
ment. Let 

x,{j) =ln[/ex(i)] (6) 

be the yearly maximum value of x in the jth year. From the measured 
data of Xi(j), j = 1, 2, 3, •••, M, we can obtain an approximate 
distribution of X\. The parameters a and U can be estimated by fitting 
the theoretical P(xi ^ X) to the measured data. Gumbel 6 has shown 
that a least-square fit of P(x\ ^ X) to the data leads to the following 
formulas for calculating a and U: 



a = — 



(7) 
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Fig. 4 — Comparison of 20-year rain-rate distribution calculated by extreme 
statistics method (dashed line) with 20-year data (solid line) for Pittsburgh, Pennsyl- 
vania, and Chicago, Illinois. 
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Fig. 5 — Comparison of 20-year rain-rate distribution calculated by extreme 
statistics method (dashed line) with 20-year data (solid line) for Harrisburg, Pennsyl- 
vania, and New Orleans, Louisiana. 
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Fig. 6 — Comparison of 20-year rain-rate distribution calculated by extreme 
statistics method (dashed line) with 20-year data (solid line) for Lynchburg, Virginia. 



and 



»-h&™*' 



= Vz 2 - z . 



(14) 



(15) 



Thus, we have all the necessary formulas for calculating P(xk ^ x). 
To obtain the rain-rate distribution, we substitute (2) into (3) to yield 



where 



4-1 p-Ny 

P(R k fc r) « l - «-«-» £ ^r , 

7V = o -AM 



(16) 



y = a[(ln r) - UJ (17) 

Therefore, the time that R k will exceed the threshold r, on long-term 



average, is 



T{R k >r) = r X P(R k ^ r) 



k-\ p -Nv 

-- r x J 1 - r** E - m - 



(18) 



where r = 5 minutes is the rain-gauge integration time. Furthermore, 
in any given year, .ft* and R t will never occur in the same 5-minute 



RAIN RATE 1117 






10 J 


- 






LA GUARDIA AIRPORT, NEW YORK 


10 2 
10 1 








\\ 


10° 










10-' 




I I 


I I 


iii 



20 40 60 



80 100 120 140 160 180 200 220 240 260 280 300 
5-MINUTE RAIN RATE IN MM/HR 



Fig. 7 — Comparison of 49-year (1903-1951) distribution of 5-minute rain rates 
calculated by extreme statistics method with measured 20-year (1953-1972) data at 
La Guardia Airport, New York, New York. 

interval if k ^ I. This means T{R k ^ r) for various order A; can be 
summed to yield an approximation to the original rain-rate distribution 
in the extremal region; i.e., 

T(R ^ r) ~ £ T(R k ^ r) (19) 

for high rain rates. 

The only input required for the calculation is the yearly maximum 
5-minute rain rates, 

RiU),3 = 1,2,3, >-,M, 
which can be obtained from National Climatic Center publications. 4 

III. COMPARISON OF CALCULATED AND MEASURED RESULTS 

Figures 1 through 6 display the comparison of the distributions of 
high rain rates calculated via statistics of extremes from the data for 
eight locations. The number of years M is 20 (from 1953 to 1972). In 
these figures, the solid lines represent the data obtained by the method 
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described in Refs. 2 and 3, whereas the dashed lines are the distribution 
calculated by the theory of extremes. 

In Fig. 3, the value of S in eq. (19) is varied from 1 to 12. For S = I, 
T(Ri ^ r) is the distribution of yearly maximum 5-minute rain rate 
and is approximately equal to T(R ^ r) only in the extremal region 
(beyond 160 mm/hr). As r decreases, T{R\ ^ r) deviates significantly 
from T(R ^ r), limiting at the 5-minute-per-year level as r approaches 
zero. The basis for this saturation is that there is only one yearly 
maximum 5-minute rain rate (with 5-minute duration by definition) 
in any given year. It is obvious that the yearly maximum 5-minute 
rain rate can exceed any threshold by no more than 5 minutes per 
year. Similarly, for S = 2, T(Rr ^ t>) plus T(R 2 ^ r) is limited to a 
10-minute-per-year level as r approaches zero. However, Fig. 3 shows 
that the applicable range of approximation (19) increases rapidly 
with S. 

For engineering terrestrial radio paths, we are interested in the 
range of rain-rate distributions below 50 minutes per year, because a 
single radio hop outage exceeding 50 minutes per year is considered 
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Fig. 8 — Comparison of 50-year (1897, 1903-1951) distribution of 5-minute rain 
rates calculated by extreme statistics method with measured 20-year (1953-1972) 
data for Pittsburgh, Pennsylvania. 



RAIN RATE 1119 




20 40 60 



80 100 120 140 160 180 200 220 240 260 280 300 

5-MINUTE RAIN RATE IN MM/HR 



Fig. 9— Comparison of 36-year (1903-1912, 1926-1951) distribution of 5-minute 
rain rates calculated by extreme statistics method with measured 20-year (1953-1972) 
data at Midway Airport, Chicago, Illinois. 

undesirable. Figure 3 shows that using S = 12 in eq. (19) includes 
the range of interest to terrestrial radio-path engineering. 

Figures 1 to 6 show similar close agreement between the calculated 
results and the data. 

The advantage of using the extreme statistics method is its sim- 
plicity; i.e., it is much easier to obtain the extreme-statistics results 
from the set of 20 maximum yearly rain rates than to obtain the other 
displayed results which require analyzing all heavy rainfalls in each 
year. 

IV. FIFTY-YEAR DISTRIBUTIONS OBTAINED FROM RAINFALL 
INTENSITY-DURATION-FREQUENCY CURVES 

Rainfall intensity-duration-frequency curves published by the 
Weather Bureau 16 were obtained from approximately 50 years (1900- 
1950) of rain-rate data processed in accordance with the statistical 
theory for distribution of yearly maximum rain rate by the Gumbel 
method. 6 In other words, these curves represent P(Ri ^ r), where Ri 
is the yearly maximum rain rate. The return period Q{R\ ^ r), which 
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is the reciprocal of P(Ri ^ r), is labeled on each curve; i.e., 

Q(Ri *r) = — ^ years. (20) 

These curves cover the range-of-return period from 2 to 100 years and 
rain-gauge integration time t (i.e., the duration) from 5 minutes to 
24 hours. We considered only the 5-minute rain-rate statistics in 
this report. 

Since these curves represent P(R\ ^ r), the two parameters a and 
U in eq. (17) can be estimated by fitting the theoretical P(Ri ^ r) 
to two points on the intensity-duration-frequency curves. Once a and 
U are obtained, we can calculate a distribution of high rain rates by 
eqs. (18) and (19). However, one adjustment is needed in this process. 
Gumbel 6 has shown that as the number of years M approaches infinity, 
Z and <r t in eqs. (13) and (15) approach the following asymptotic values 

lim Z = 7 = Euler's Constant 
*-•• ~ 0.5772 (21) 

lim a. = ~ (22) 

The corresponding asymptotic forms for a and U are 



and 



lim a = a a = -£- , (23) 



lim V - V„ = xi - ^-- (24) 



The rainfall intensity-duration-frequency curves were obtained using 
the asymptotic values, a„ and U„, even though the number of years 
M are 50 or less. Such approximations introduce slight errors and can 
be corrected by the following relationships among a m , [/„, a, and U: 

a = a^-Oz (25) 



*-«■+={»-!■*] 



(26) 



To relate theoretical P(Ri ^ r) with the intensity-duration-frequency 
data, we combine eq. (16) and (20) to give 

Q(Ri fc r) - 1 _ *_,-,, , (27) 

where 

2/ = a.Cdn r) - 17.]. (28) 
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The two parameters, a M and [/«,, can be determined by equating (27) 
to two sets of data, say (r„, Q a ) and (r b , Q b ), read from the intensity- 
duration-frequency curves. It is easily shown that the relationships are 

A, -A t (29) 

In r a — In n 
and 

. _ A a In r b — A b In r a f9ny 

U « ~ A a -A b ' {dU) 



where 



and 

A b = —In 



„„ = _ ln [ ln _%_], (31 ) 

[ In o^r]- (32 > 



By substituting a M and £/„ into (25) and (26), we obtain a and U for 
the 50-year data. Substituting a and U into (18) and (19) gives the 
50-year distribution of 5-minute rain rates. 

The time bases (i.e., M) in the intensity-duration-frequency curves 
are mostly 50 years or less. However, several locations have time bases 
much shorter than 50 years. For example, the time bases are 18, 16, 
and 17 years for Mt. Tamalpais, California; Tonopah, Nevada; and 
Yakima, Washington, respectively. Due to this limitation, we have 
chosen Q a = 2 years and Q b = 10 years for calculations of a and U. 
Since Q a and Q b are fixed, we then need only three numbers: M, r a , 
and r b for each location, read from the intensity-duration-frequency 
curves to calculate the rain-rate distribution. 

For example, for New York City, the three numbers are : 

M = 49 years (1903-1951) 

r a = 4.4 inches/hr = 111.8 mm/hr 

r b = 6.5 inches/hr = 165 mm/hr. 

Substituting r a and r b into eqs. (29) and (30) yields 

<*„ = 4.828 
U„ = 4.64. 

Substituting M into eqs. (12) through (15) and a M and £/„ into (25) 

and (26) gives 

a = 4.363 

U = 4.63. 

The 49-year (1903-1951) distribution of 5-minute rain rate calculated 
from this a, U pair for New York City is very close to the 20-year 
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(1953-1972) data, as shown in Fig. 7. Figure 8 shows similar close 
agreement between 50-year and 20-year distributions in Pittsburgh, 
Pennsylvania. 

However, Fig. 9 shows an appreciable difference between 36-year 
and 20-year distributions for Chicago, Illinois. This appreciable 
difference, the irregular shape of the 20-year distribution in Fig. 9 and 
the instability noted in Figs. 1 and 2 indicate that a 20-year time base 
with a single rain-gauge measurement may not be sufficient to guaran- 
tee a stable distribution for extremely high rain rates. If more stable 
results are required, the 20- and 50-year data may be combined to give 
a 70-year distribution. 

The calculated curves (dashed lines) in Figs. 1 through 6 are based 
on 20-year data from 1953 to 1972 in Ref. 4, whereas the calculated 
curves (dash-dot lines) in Figs. 7 through 9 are based on approximately 
50 years of data in Ref. 16. 

V. CONCLUSION 

A method has been described for calculating the distribution of high 
rain rates by applying the statistical theory of extremes to the available 
yearly maximum 5-minute rain-rate data. Figures 1 through 9 show 
that the calculated distributions agree closely with the data in the 
heavy-rain region of interest to radio-path engineering. The virtue of 
this method is that only yearly maximum rain-rate data are required 
to generate satisfactory results for radio-path engineering. The rainfall 
intensity-duration-frequency curves 16 provide approximately 50 years 
(1900-1950) of such data for 203 locations in the United States. 
Furthermore, a new publication on "Maximum Short-Duration Pre- 
cipitation," 17 for approximately 300 U. S. locations issued annually by 
the National Climatic Center since 1973, provides additional yearly 
maximum rain-rate data. Therefore, long-term (^50 years) distri- 
butions of high rain rates for 203 U. S. locations can easily be obtained 
by this method. 
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